On Defect-Mediated Transitions in Bosonic Planar Lattices 
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Abstract 

We discuss the finite-temperature properties of Bose-Einstein condensates 
loaded on a 2D optical lattice. In an experimentally attainable range of 
parameters the system is described by the XY model, which undergoes a 
Berezinskii-Kosterlitz-Thouless (BKT) transition driven by the vortex pair 
unbinding. The interference pattern of the expanding condensates provides 
the experimental signature of the BKT transition: near the critical tempera- 
ture, the k = component of the momentum distribution sharply decreases. 



Bose-Einstein condensates (BECs) are nowadays routinely stored in ID [1-5], 2D [6] and 
3D [7] optical lattices. In this paper we review and further discuss the finite-temperature 
properties of Bose-Einstein condensates loaded on a 2D optical lattice [8]: there we showed 
that at a critical temperature lower than the temperature Tbec at which condensation 
in each well occurs, 2D lattices of BECs may undergo a phase transition to a superfluid 
regime where the phases of the single-well condensates are coherently aligned allowing for 
the observation of a Berezinskii-Kosterlitz-Thouless (BKT) transition. 

The BKT transition is the paradigmatic example of a defect-mediated transition [9] which 
is exhibited by the two-dimensional XY model [10,11], describing 2-components spins on 
a two-dimensional lattice. In the low-temperature phase, characterized by the presence of 
bound vortex-antivortex pairs, the spatial correlations exhibit a power-law decay; above 
a critical temperature T BKT , the decay is exponential and there is a proliferation of free 
vortices. The BKT transition has been observed in superconducting Josephson arrays [12] 
and its predictions are well verified in the measurements of the superfluid density in 4 He 
films [13]. In finite magnetic systems with planar symmetry [14], the BKT transition point 
is signaled by the drop of a suitably defined magnetization [15]. We shall discuss in the 
sequel the analogous of this magnetization for atomic systems. 

For 2D optical lattices, when the polarization vectors of the two standing wave laser fields 
are orthogonal, the resulting periodic potential for the atoms is V opt = Vo[sin 2 (kx)+sin 2 iky)] 
where k = 2n/\ is the wavevector of the laser beams. The potential maximum of a single 
standing wave Vq = sEr may be controlled by changing the intensity of the optical lattice 
beams, and it is conveniently measured in units of the recoil energy E R = h 2 k 2 /2m (m is the 
atomic mass), while, typically, s can vary from up to 30. Gravity is assumed to act along 
the z-axis. Usually, superimposed to the optical potential there is an harmonic magnetic 
potential V m = {m/2)[uj 2 {x 2 + y 2 ) + oj 2 z 2 }, where oj z {uj r ) is the axial (radial) trapping 
frequency. The minima of the lattice are located at the points j = (ji, J2) • § with ji and 
j 2 integers, and the potential around the minima is V opt ~ ^uj 2 \^x — Aji/2) 2 + (y — Aj'2/2) 2 ] 
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With Ld r = ^2V k 2 /m. When uj r ^> uj r ,uj z , the system realizes a square array of tubes, i.e. 
an array of harmonic traps elongated along the z-axis [6]. In the following we analyze only 
the situation in which the axial degrees of freedom are frozen out, which is realizable if u z 
is sufficiently large. We also assume that the harmonic oscillator length of the magnetic 
potential \Jh/ muj r is enough larger than the size L (in lattice units) of the optical lattice 
in order to reduce density inhomogeneity effects and to allow for a safe control of finite size 
effects. 

When all the relevant energy scales are small compared to the excitation energies one can 
expand the field operator [16] as ^(f, t) = Y^jipj(t)&j(r) with $j(r) the Wannier wavefunc- 
tion localized in the j-th well (normalized to 1) and = tpt^ the bosonic number operator. 

J J J 

Substituting the expansion in the full quantum Hamiltonian describing the bosonic system, 
leads to the Bose-Hubbard model (BHM) [17,16] 

H = -K g + /i.c.) + E 

<?J> 3 

where E<i j> denotes a sum over nearest neighbours, U = (AiiTi 2 a/m) J dr&i (a is the s-wave 
scattering length) and K ~ - / dr [£v$ ? • V$j + ^V ext ^>.]. 

Upon defining J 2KN (where A^ is the average number of atoms per site), when 

A^ > 1 and J/N$ < U, the BHM reduces to H = H X y - f Ejft, which describes the 

i 

so-called quantum phase model (see e.g. [18,19]): 9j is the phase of the j-th condensate and 

Hxv = -J E cos(%-%) (2) 

<?J> 

stands for the Hamiltonian of the classical XY model. When J ^> U and at temperatures 
T ^> U/Icb, the quantum phase Hamiltonian may be well approximated by Eq.(2); therefore, 
the pertinent partition function describing the thermodynamic behaviour of the BECs stored 
in an optical lattice can be computed in these limits with the classical XY model. 

We remark that the BHM - for U = - describes harmonic oscillators and thus cannot 
sustain any BKT transition. However, in the limits above specified (iVo ^> 1 and J/Nq -C 
U J) the Bose-Hubbard Hamiltonian reduces in the XY model (2), which displays 
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the BKT transition. We observe that in the chosen range of parameters, the condition 
UNq/ J = UN /2K 3> 1, satisfies the finite-temperature stability criterion recently derived 
in [20]. In superconducting networks N is the average number of excess Cooper pairs and 
thus, the requirements N 3> 1 and J/Nq <C U are always satisfied; at variance, in bosonic 
lattices N varies usually between ~ 1 and ~ 10 3 and the validity of the mapping in the 
quantum phase Hamiltonian is not always guaranteed. 

A simple estimate of the coefficients J and U may be obtained by approximating the 
Wannier functions with gaussians, whose widths are determined variationally [8]. For a 2D 
lattice with V between 20 and 25E R (and N m 170 as in [6]), the conditions J > U > J/N$ 
are rather well satisfied and the BKT critical temperature, T BKT ~ J/k B , is between 10 and 
30nK. 

We observe that, using a coarse-graining approach to determine the finite-temperature 
phase-boundary line of the Bose- Hubbard model (1) [21], one gets - in 2D for N 3> 1 and 
J/U ^> 1 - a critical temperature T BKT m J/k B . Indeed, within a coarse-graining approach, 
the phase-boundary line of the Bose-Hubbard model at a temperature T is given by [21] 

1 = Kd [ drG(r) (3) 
Jo 

where d is the lattice dimension and (3 = l/k B T. In Eq.(3) G(r) is the correlation function, 
which, for N » 1, is given by 



G(t) » iV ^^ J e -TU(l-r/p)/2 (4) 



^3 \nN ,q 

with & 3 (z,q) = l + 2J2n=i C0S ( 2nz )-Q n2 - Since d = 2, and J = 2KN , from Eq.(3) it follows 
1 = JG(uj = 0)/N where G(u = 0) = $g drG(r); one has then 



G(u = 0) _ E m exp (-pUft/2) ■ l - c %^ m 



(5) 



No E m exp(-^02/2) 

where <p m = N —m. Eq.(5) yields, for U <ti k B T, G(u = 0)/N pz (3, thus implying (3 J m 1, 
i.e. k B T BKT pz J, in fair agreement with XY estimates. 



We notice that the XY Hamiltonian, in the limits N ^> 1 and J 3> U 3> J/Nq, 
can be retrieved extending the Gross-Pitaevskii (GP) Hamiltonian to finite temperature 
T U/ks- In fact, replacing the tight-binding ansatz \I/(r, t) = Y^jipj(t)^j(r) (where ipj 
are classical fields and not operators) in the GP Hamiltonian one gets the lattice Hamiltonian 
H = -KY, < X3 > hl^'^j + c - c -) + ( u / 2 ) T,jNj(Nj - 1), which is the classical version of the 
Bose-Hubbard model (Nj =| ipj | 2 ). Writing Nj — N + SNj, one may neglect the quadratic 
terms [i.e. (5Nj) 2 } in the hopping part of the Hamiltonian —KY^^j^ipiipj + c.c), which 
then reduces to H X y (2). 

Accurate Monte Carlo simulations yield - for the XY model - the BKT critical tempera- 
ture T BKT = 0.898 J/k B [22]. When U <C J, a BKT transition still occurs at a slightly lower 
critical temperature T B kt(U). To evaluate more accurately the effects of the interaction 
term on the BKT critical temperature one may use a semiclassical approximation to eval- 
uate a renormalized Josephson energy [23,24]: the renormalization group equations yields 
then the following equation for /C = (3T BKT (U) 

1 1 



2-tt/C 



where X u = U/nJ and g()C,X u ) = E^°=i[A u /C 2 /2 - (n 2 /2)ln(l + X u JC 2 /n 2 )}. One may see 
[25] that for U / J < 36/V, Eq.(6) has an unique solution; at variance, for U / J > 36/7T Eq.(6) 
does not have any solution. The value U/J = 36 /n is in reasonable agreement with the 
known mean-field prediction 2zJ/U xs 1 for the T = phase transition (see, e.g., Eq. (30) 
of Ref. [26] for A^ ^> 1). In the inset of Fig.l we plot the value of T B kt(U)/T B kt (where 
Tbkt is the BKT critical temperature of the XY model with U = 0) as a function of U/ J 
from the numerical solution of Eq.(6). 

The emerging physical picture is the following: There are two relevant temperatures 
for the system, the temperature T BEC , at which in each well there is a condensate, and 
the temperature T B kt at which the condensates phases start to coherently point in the 
same direction. Of course, in order to have well defined condensates phases one should 
have T BK T < Tbec- The critical temperature T B ec ~ 0.94hN^ 3 u/kB where oo = (d^c^) 1 / 3 . 



With the numerical values given in [6], T B ec turns out to be > 500nK for s > 20. When T < 
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T BEC , the atoms in the well j of the 2D optical lattice may be described by a macroscopic 
wavefunction ipj. Furthermore, when the fluctuations around the average number of atoms 
per site No 1 are strongly suppressed, one may put, apart from the factor constant 
across the array, ipj = e 9 ^ . The temperature Tbkt is of order of J / k B ' with the experimental 
parameters of [6] and V between 20 and 25E R , one has that T BKT is between 10 and 30nK, 
which is sensibly smaller than the condensation temperature Tbec of the single well. A 
similar picture describes also 2D arrays of superconducting Josephson junctions [18,19]: they 
exhibit a temperature Tbcs at which the metallic grains placed on the sites of the array 
become (separately) superconducting and the Cooper pairs may be described by macroscopic 
wavefunctions. At a temperature T B kt < Tbcs, the array undergo a BKT transition and 
the system - as a whole - becomes superconducting. 

The experimental signature of the BKT transition in bosonic planar lattices is obtained 
by measuring, as a function of the temperature, the central peak of the interference pattern 
obtained after turning off the confining potentials [8]. In fact, the peak of the momentum 

— * 

distribution at k = is the direct analog of the magnetization of a finite size 2D XY magnet. 
For the XY magnets, the spins can be written as Sj = (cos 9j, sin 9j) and the magnetization 
is defined as M = (1/iV) ■ ( | Y^jSj \ ) where (• • ■) stands for the thermal average. A spin- 
wave analysis at low temperatures yields M = (2N)~ ksT/8wJ [14,15]. With discrete BECs 
at T = 0, all the phases 9j are equal at the equilibrium and the lattice Fourier transform of 
V'j) i'k = Jr ^1^] e ~ lk " ') exhibits a peak at k — (k is in the first Brillouin zone of the 2D 
square lattice) and the magnetization is: 

M=(|^o|). (7) 

The intuitive picture of the BKT transition is then the following: at T = 0, all the spins 
point in the same direction. Increasing the temperature, bound vortex pairs are thermally 
induced. In Fig. 2, left inset, a single free vortex is depicted: we plot a configuration of the 
phases such that 9j equals the polar angle of the site j with respect to the core vortex (notice 
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that with 9 j equals the polar angle plus 7r/2 one would have the usual plot of a vortex). As 
one can see, a single vortex modifies the phase distribution also much far from its core and the 
square modulus of its lattice Fourier transform | ip£ | 2 has a minimum at k — 0. At variance, 
a vortex-antivortex pair [see Fig.2, right inset] modifies the phase distribution only near the 
center of the pair (in this sense is a defect) and its lattice Fourier transform has a maximum 
at k — 0. Increasing the temperature, vortices are thermally induced. For T < T BKT only 
bound vortex pairs are present, and in average the spins continue to point in the same 
direction. When the condensates expand, a large peak (i.e., a magnetization) is observed 
in the central k = momentum component. Rising further the temperature, due to the 
increasing number of vortex pairs, the central peak density decreases from the T = value. 
For T pa T BKT the pairs starts to unbind and free vortices begin to appear, determining a 
sharp decrease around T BKT of the magnetization. At high temperatures, only free vortices 
are present, leading to a randomization of the phases and to a vanishing magnetization. In 
Fig.l we plot the intensity of the central peak of the momentum distribution (normalized 
to the value at T = 0) in a 2D lattice as a function of the reduced temperature k B T/J, 
evidencing the sharp decrease of the magnetization around the critical temperature [8] . 

— * 

At t = the amplitude of the k = peak of the momentum distribution is simply given 
by the thermal average of ipo- By measuring the k — peak (i.e., ( | ip | 2 )) at different 
temperatures, one obtains the results plotted in Fig.l. The figure has been obtained using 
a Monte Carlo simulation of the XY magnet for a 40 x 40 array: we find k B T c ~ 1.07 J. 
In Fig.l we also plot the low-temperature spin wave prediction [15] (solid line). At times 
different from t — 0, the density profiles are well reproduced by the free expansion of the ideal 
gas. One obtains 4>(p,t) = x(Pz)e- i[iPz+m9t)S - p " ]/6m29h ^(p x ,p y ,t) [where X (Px) oc e~ a * p *' 2h2 
and similarly for xiPy) an d x{Pz)}, i-e. a uniformly accelerating motion along z and a 
free motion in the plan x — y, with <p(p x ,Py,t) = x(Px)x(Py)' l Pk e ~ l( ' p * +p ^ t ^ 2hm giving the 
central and lateral peaks of the momentum distribution as a function of time for different 
temperatures. 

An intense experimental work is now focusing on the Bose-Einstein condensation in two 
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dimensions: at present a crossover to two-dimensional behaviour has been observed for Na 
[27] and Cs atoms [28]. Our analysis relies on two basic assumptions; namely, the validity 
of the tight-binding approximation for the Bose-Hubbard Hamiltonian and the requirement 
that the condensate in the optical lattice may be regarded as planar [29]. It is easy to see 
that the first assumption is satisfied if, at T = 0, V 3> // (where \x is the chemical potential), 
and, at finite temperature, huj r > k B T. The second assumption is much more restrictive 
since it requires freezing the transverse excitations; for this to happen one should require a 
condition on the transverse trapping frequency u z . Namely, one should have that, at T = 0, 
Tiuj z ^ 8K and that, at finite temperature, Tluj z > k B T (since uj z <C u r , the latter condition 
also implies that hCu r > k B T). In [6] it is V = s ■ k B ■ 0.15/ifT and for s > 20 the tight- 
binding conditions are satisfied since 10Hz > 8K/2nh; furthermore, if u z = 2n ■ 1kHz, one 
may safely regard our finite temperature analysis to be valid at least up to T ~ 50nK. We 
notice that the experimental signature for the BKT transition for a continuous (i.e., without 
optical lattice) weakly interacting 2D Bose gas [30] is also given by the central peak of the 
atomic density of the expanding condensates. 

Our study evidences the possibility that Bose-Einstein condensates loaded on a 2D op- 
tical lattice may exhibit - at finite temperature - a new coherent behaviour in which all the 
phases of the condensates located in each well of the lattice point in the same direction. The 
finite temperature transition, which is due to the thermal atoms in each well, is mediated by 
vortex defects and may be experimentally detectable by looking at the interference patterns 
of the expanding condensates. Our analysis strengthens - and extends at finite temperature 
- the striking and deep analogy of bosonic planar systems with 2D Josephson junction arrays 

[!]■ 
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FIGURES 
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k B T/J 

FIG. 1. Intensity of the central peak of the momentum distribution (normalized to the value 
at T = 0) as a function of the reduced temperature t = ksT/J; empty circles: Monte Carlo simu- 
lations; solid line: low-temperature spin wave prediction [15]. Inset: the BKT critical temperature 
Tbkt(U) (in units of the XY critical temperature) as a function of U/J from Eq.(6). 



FIG. 2. 
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